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Nearest-Neighbour Contributions 
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Superconducting phase diagram of the extended Hubbard model supplemented with interaction 
and hopping terms exceeding nearest neighbour distance in range is analysed systematically at differ- 
ent band-filling and temperature values in a mean-field approximation. The obtained results clearly 
underline the importance of next-nearest neighbour terms in developing the main superconducting 
properties of the model system. In particular, the emergence of superconducting phases of different 
■ symmetry at a given point of the phase diagram, critical temperatures T c , zero temperature gap 

amplitude values Ao, Ao/T c ratios, doping and temperature dependences are all strongly influenced 
by next-nearest neighbour contributions. 



■ I- INTRODUCTION 

H 

The subject of high critical temperature superconductors starting from an early proposal of Anderson (Anderson, 
1987) polarizes significant forces around the Hubbard model, its extensions and derivatives (Dagotto, 1994). On 
C$ ■ this line, taking the fact into account that high temperature superconductivity generally emerges in systems with 
a pronounced planar structure (Smith, Manthiram, Zhou, Goodeneough and Market, 1991), the two-dimensional 
square lattice version of the Hubbard model in its extended version is intensively used as a starting point for describing 
superconducting properties of this type of materials (Micnas, Ranninger and Robaszkiewicz, 1990). As a consequence, 
within the last few years the Hubbard Hamiltonian containing nearest neighbour (AW) interaction and hopping terms, 
has been extensively studied (de Boer, de Chatel, Franse and de Visser, 1995) and the emergence of superconductivity 
has been rigorously demonstrated in some conditions (de Boer, Korepin and Schadschneider, 1995) in these systems. 
In the last period, however, a great number of results directed our attention to the importance of extending the 
interaction and hopping range in such type of an analysis. The stimulating observations can be summarized as 
follows. 

£ — ' Concerning the kinetic energy term, Koma and Tasaki starting from the approximation-free study of two-point 
correlation functions at finite temperatures in low dimensional extended Hubbard models containing more than AW 
00 ■ hopping terms (Koma and Tasaki, 1992) sharply demonstrated the fact that electron hopping with higher than AW 
terms plays a fundamental role in various condensation phenomena in itinerant-electron systems. The importance of 
hopping-range in building up the phase diagram was also emphasized in an exact manner for finite- (Tasaki, 1995) 
and infinite- U (Verges, Guinea, Galan, van Dongen, Chiappe and Louis, 1994) Hubbard models. Furthermore, as 
it was pointed out by Veilleux et al. (Vcillcux, Dare, Chen, Vilk and Tremblay, 1995) within a model containing 
only NN hopping the Fermi surface topology and the filling dependence of both the Hall coefficient and the uniform 
magnetic susceptibility are qualitatively wrong, a disagreement that cannot be removed perturbatively. On the other 

■ hand, incorporating also next-nearest neighbour (NNN) hopping the band structure becomes more realistic and all the 
above mentioned physical quantities, as well as the position of neutron scattering intensity maxima, have the correct 
qualitative behaviour (Lavagna and Stemmann, 1994, Littlewood, Zaanen, Aeppli and Monien, 1994). Following 
this line of reasoning it is important to mention that the relative influence of density of state (DOS) effects - like 
van Hove singularities (Blumberg, Stojkovic and Klein, 1995) - and Fermi surface topology effects - like nesting 
(Gulacsi, Bishop and Gulacsi, 1995) - on pairing are extremely important in describing superconducting properties 
of high-Tc materials. While both van Hove singularity and nesting occur simultaneously at half filling in the usual 
AW model, there is no nesting when NNN hopping is present (Veilleux et al. and references cited therein) and the 
influence of saddle point effects (Abrikosov, 1995, Abrikosov, 1994) can be taken more accurately into consideration. 
Furthermore, the presence of hopping exceeding AW terms in range was observed, in particular, by thermopower 
measurements in related systems (Ponnambalam and Varadaraju, 1995) and was found to be consistent with angle- 
resolved photoemission data (Fehrenbacher and Norman, 1995) and gap symmetry analysis (O'Donovan and Carbotte, 
1995) too. We also would like to mention that the ratio of AW to NNN hopping amplitudes is clearly known for a 
broad spectrum of high-Tc materials (Brenig, 1995). 

Concerning the interaction terms, it is known from exact studies of phase diagrams related to extended Hubbard 
models in higher than one dimension (Strack and Vollhardt, 1994) that contributions exceeding AW distances in 
range can significantly influence the emergence of condensed phases. Up to this moment such type of studies have 
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been carried out on magnetic phases, however, the relevance of the obtained results to superconductivity is quite 
straightforward (Strack and Vollhardt, 1995). Similarly, the non-negligible character of NNN interaction contributions 
was clearly pointed out experimentally. For example, the importance of NNN type off-site interactions was revealed in 
the interpretation of Auger core- valence line shapes (Verdozzi and Cini, 1995), which allows the direct observation and 
measurement of interparticle interaction strength and its distance dependence in valence bands of solids (Verdozzi, 

1993) . Furthermore, King et al. based on angle resolved photoemission data (King et at, 1993), clearly underline that 
model Hamiltonians containing interaction terms only up to the NN ones are not sufficient for the proper description 
of the electronic structure near the Fermi energy (Ep) in NCCO type materials. They also conlude that NNN 
interactions are very important and models including only NN interactions in Cu0 2 planes are unable to explain 
related clear experimental facts. The interpretation of two-magnon Raman scattering in high-T c materials (specially 
the A\ g and B 2g finite intensity peaks) also requires second-nearest neighbour or even longer-range contributions to 
be incorporated at the level of the interaction terms (Brenig). Moreover, the importance of NNN terms in developing 
basic properties of extended Hubbard models has been claimed (Di Stasio and Zotos, 1995), their effect on screening 
processes and as well as on phase diagram configurations has been studied (van den Brink, Meinders, Lorenzana, 
Eder and Sawatzky, 1995), the relevance to YBCO materials has been emphasized for various conditions (Grigelionis, 
Tornau and Rosengren, 1996) and their importance in loss of antiferromagnetic order has been discussed (Kampf, 

1994) . 

Stimulated by these findings we present an extensive study of superconducting phase diagram of the extended 
Hubbard model containing more than nearest neighbour contributions, in order to provide a systematic image about 
the importance of the longer range contributions. The study is given in two dimensions, based on a traditional 
Hartree-Fock type decoupling procedure in terms of the Green's function description allowing, however, arbitrary gap 
symmetry. The DOS is treated exactly, and the stable superconducting state is chosen on the basis of a free-energy 
type analysis. Such type of a systematic investigation, incorporating the effects of NNN terms too, has not been 
presented so far in literature. 

Our results emphasize the importance of interaction and hopping terms exceeding nearest neighbouring distance in 
range in building up superconducting characteristics of the system described by an extended Hubbard type model. 
In particular, the emergence of superconducting phases of different symmetry at a given point of the phase diagram, 
critical temperatures, zero temperature gap amplitude values, Aq/T c ratios, doping and temperature dependences are 
strongly influenced by the presence of NNN terms. Based on the results we conclude that the effect of NNN terms is 
not marginal but essential in the description of basic superconducting properties. 

Th e pap er is organized as follows. Section [II] presents the model, Sec. EH contain s besi des the numerical procedure 
(Sec. [II A ) the obtained results grouped into three paragraphs as follows: Sec. [II B analyzes the AW case, Sec. 
Ill C outlines the effect of hopping terms exceeding AW distances in range and Sec. Ill D presents the effects of NNN 
interaction contributions, the study being given at various temperature and doping values. Finally, the conclusions 



and summary close the presentation in Sec. IV 



II. THE MODEL AND ITS DESCRIPTION 



In this paper the following extended Hubbard Hamiltonian is considered on a square lattice: 

ff = -EE* y 4 c 3a + ^ E E w' n ™ n J*' > (!) 

i,j a i,j a.a' 

where the fermion operators c\ a and Ci a create and annihilate, respectively, electrons with spin a in the single 
tight-binding orbital associated with site i, and riic is the particle number operator. The parameters ti = Uj with 
I = 1, 2, ... ,5 are the real space hopping matrix elements between sites i and j with a diagonal part to = ta. The 
coupling constants are independent parameters denoted below as Uq = U^~ a , U\ = UfJ +1 and U% = U[J +2 

on-site, nearest neighbour and next-nearest neighbour contributions, respectively. All higher order interaction terms 
(for j >i + 2) will be neglected for the present study, that is, Ui = U[J >i+2 — 0. Furthermore, any of Uq, U\ and U 2 
can either be positive or negative with no further restriction. This general choice of interaction constants allows us a 
systematic exploration of the full superconducting phase diagram of the model under study. 

The gap equation obtained within a Hartree-Fock decoupling procedure (Abrikosov, Gor'kov and Dzyaloshinskii, 
1965) is given by 

A(fc) = - U ($- k) m F(q) , (2) 
<7 
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and the chemical potential // is controlled via 



1 - S = £ ( 1 + m F(k) ) . (3) 
% 

Here (5 = 1 — 2n is the band filling factor with positive or negative values corresponding to hole or electron doping, 
respectively. Furthermore, n is the particle number density for a given spin index. In the above equations 

= tanh(l^(fc)) 
2E(k) 



E{k) = yj W{k) + A^k) , (5) 
fl(k) = fi - e(k) - nU, (6) 

with U — Uo + 8Ui + 8U2, and (3 stands for the reciprocal temperature. Taking the lattice spacing to be unity, the 
explicit expression of the inter-particle coupling in k-space for a square lattice is given by 

U(k) = Uq + 2U\ (cosk x + cosfc. y ) + AU2 cosfc^cosfcy . (7) 
One can easily verify that in Eq. (||) 

U{q-k)= 9iMk)%(q), (8) 

i G TZ 

where TZ collects all the irreducible representations of the point group (Cn v ) characterizing the square lattice, ^i(k) 
represents the basis functions of the ith irreducible representation and gi denotes the effective coupling constants (see 
Table Q). As a consequence, the k-dependent order parameter can be given in terms of these basis functions, that is, 

A(fc) = A * • ( 9 ) 

i e K 

In Eq. (||) i represents both singlet and triplet pairing. Based on Eqs. (||) and (Q) the gap amplitudes can be written 
as 

A 4 = -g i ]T Aj F ® Mk) (10) 

for each irreducible representation i. In case of representations B\ and B2 the gap equation is scalar, while for 
representation A\ we have to solve a coupled system consisting of three equations due to mixing of the corresponding 
three basis functions. A similar situation emerges for triplet pairing where two equations are coupled together since 
the corresponding irreducible representation of the states are two-dimensional. 



TABLE I. The basis functions 9i(k) and their notations, representations and the corresponding effective coupling constant 
gi entering the pair potential. 



TZ 

Ai 
Ax 
Si 
B 2 



notation 



1 

i (cos k x + cos 
cos k x cos k y 



ky) 



y 

|(cos k x — cos k y ) 
sin k x sin k y 

sin k y 
sin k y 



( sin k x - 
I sin k x ■ 
k x cos k y + sin k y cos k x 
kxCos ky sin ky cos kx 



s 

s* 

$ x y 
x 2 —i 

P 
p 

P2 
P2 



U 
4Ui 
W 2 

4(7i 
8U 2 
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The dispersion relation is given in terms of the real space hopping matrix elements as 



e(fc) = t + ^2 U exp(ikR lii+ i) , (11) 
i 

where Ri t i+i denotes lattice vector pointing from site i to its I th neighbour (l> 1). 

At a given temperature all the possible ordered phases can be found as the nontrivial solutions of Eqs. (|2|) and 
(||). After the self-consistent solution, the stability of ordered phases of different symmetry was carefully investigated. 
This has been performed on the basis of a comparative free-energy (or at T — internal energy) analysis (for the 
procedure see eg. (Dahm, Erdmenger, Scharnberg and Rieck, 1993)). If the gap-equation presented more possible 
solutions connected to the same point of the phase diagram, the free-energy of each solution was determined and 
compared. For any solution, the free-energy relative to the normal state per particle SF, can be derived by improving 
the method of Gyorffy et al. (Gyorffy, Staunton and Stocks, 1991), the final result being given by 

SF[A(k), fj,, T] = — E(k) - \m\ } + 2 W A»(S) F(k) - i In } 

+ EE A ^ F ^ u &- *) F ® ■ ( 12 ) 

fc <? 

In this manner, altering the interaction parameters, the superconducting phase diagram of the present model has been 
constructed systematically at different temperature and band-filling values and main superconducting characteristics 
have been studied in every phase diagram domain in detail. 



III. THE OBTAINED RESULTS 



A. Numerical procedure 

To construct a phase diagram as presented above, first Eqs. @ and ( |l0| ) have been solved for the amplitudes A, 
corresponding to different representations. Since the factor F(k) occurring in these equations depends on all the Aj 
values, a self-consistent procedure has to be applied. This requirement has been fulfilled by the adoption of Broyden's 
algorithm (Numerical Recipies in Fortran, 1992), which assures the global convergence. 

To start the algorithm, first a vector composed of Aj and fx denoted by Aj 1 ' is picked up arbitrarily at the fixed 
T, n and gi values. At this choice the possibility of mixing different symmetry representations is allowed with no 
restrictions. Then using this point as a zeroth guess for the solution of Eqs. (||) and (|l0|) the nonlinear set of equations 

is solved numerically in order to get A^ 1 -*. In the knowledge of Aj the initial guess for the next iteration A^ can 

be constructed in terms of the linear combination of Aj 1 ^ and A . This linear combination provides the global 
nature of convergence (i.e. directional dependences are eliminated and the solution does not depend on the starting 
A- 1 ' value). Then, using A,- 2 - 1 as a starting point for the next iteration (to = 2), the equations are solved again 
to obtain A^ 2 \ in terms of which the starting point for the 3rd iteration A- 3 ^ can be constructed, and so on. The 
above procedure is continued until self-consistency has been reached, that is, until at the end of the mth iteration the 
difference e m = |Ai m -* — A_j m -'| is less than a small, positive and a priori fixed value, e. In the present calculation e is 

chosen to be 10~ 6 . The required convergence is fulfilled in average after to ~ 25 — 30 iterations for the chosen value 
of e. 

Because of the complexity of this analysis, only the possible superconducting phases were considered (i.e. the 
possible emergence of spin (SDW) and charge (CDW) density wave type phases was discussed, although it was not 
investigated systematically in detail). 



B. Hamiltonian with nearest neighbour terms only 



We started our study with the analysis of pure NN-case where [^2 = is satisfied and non- vanishing hopping matrix 
elements were taken into consideration only between nearest neighbouring sites, that is, ti ^ for I = 1 and <; = 
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otherwise. In what follows, all the couplings are normalized by the bandwidth, W = 8ti (for example, in the case of 
Bi2212 we have W = 1.192 eV extracted from ARPES data). The obtained results are plotted in Fig. 0. Here four 
separate plots are presented at zero (Fig. |(a) and (c)) and nonzero (Fig. 0(b) and (d)) temperatures for half filling 
(5 = 0.0, Fig. 0(a), (b)) and away from half filling (<S^0.0, Fig. 0(c), (d)). 

Concerning the notations, paramagnetic refers to the normal state of the model, d-wave and p-wave stand for order 
parameters with (cos k x cos k y ) and (sin/c^-i-sin/cy) k-dependence, respectively, while A\ denoting the ^4i-symmetry 
solution of the gap equations represents in this paragraph (in Fig. 0) a standard s-wave component for plots (a) and 
(b), and an (s,s*) admixture for plots (c) and (d). Here s* stands for the extended s-wave solution. 

In Fig. 0(a) the T — ground state phase diagram at half-filling is shown. As it can be seen, if either the on-site 
Uq or the intersite U\ interaction is attractive superconductivity emerges, where the symmetry of the stable, ordered 
phase revealed by the 8F analysis is or d x 2_ y 2 for resonable values of |J7i/W|. Furthermore, if both interaction 
parameters are attractive a competition between different symmetry pairing states is observed. It can also be seen 
that increasing the value of \Uq\, triplet pairing is no longer favored. This means particularly, that repulsive on-site 
interaction stabilizes the superconducting state of rf 3 .2_ y 2-symmctry 

Comparing these facts with the preliminary results of Micnas et al., a qualitative agreement can be observed. Here 
we should mention that taking spin and charge density wave ordering possibilities also into account, the paramagnetic 
region of Fig. 0(a) can be filled up by these phases not affecting the main superconducting regions of the plot. 
However, the occurrence of pure d x 2_ y 2 state instead of an s-d mixture (Micnas, Ranninger, Robaszkiewicz and 
Tabor, 1988, Micnas, J. Ranninger and S. Robaszkiewicz, 1989) is found in Fig. 0.a, which is in accordance with the 
results of Fehrenbacher and Norman (Fehrenbacher et al.), O'Donovan and Carbotte (O'Donovan et al.) and also 
with group-theoretical studies (Wenger and Ostlund, 1993). 
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FIG. 1. Nearest-neighbour phase diagram at zero (fcgT = 0.0) and nonzero (fesT = 0.005IF) temperatures for different 
dopings (8), using dispersion containing only nearest neighbour hopping terms. Dotted lines mean coordinate axis, solid lines 
represent boundaries between different symmetry phases. For notations A\, d, p, and paramagnetic see the text. 



With the non half-filled band case at zero temperature (Fig. 0(c) at 8 = 0.2), wc start to explore a parameter domain 
of the phase diagram which has not been investigated systematically so far in the literature. As it can be seen from the 
plot, the symmetry of the pairing state is strongly doping dependent; the A\- and p-wave solutions, with increasing 
6, axe gradually overwhelming the superconducting phase of d x 2_ y 2 symmetry, as it is predicted by Micnas et al.. On 
the other hand, moving away from half-filling an s* component of the Ai-symmetry solution sets in, suggesting that 
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the rigorous constraints on s-wave pairing deduced by Zhang for the standard Hubbard model with on-site interaction 
only (Zhang, 1990), is also applicable in the extended Hubbard model case. 

The T^O counterparts of Fig. |l](a) and (c) are Fig. [l](b) and (d), respectively, keeping the same notations for the 
labeling of superconducting phases. In Fig. |l](b) the phase diagram for the half-filled band is depicted. It can be 
observed, around |£/o/W| ss 0.15, that with increasing temperature the stable A\- and d a .2_ y 2-symmetry supercon- 
ducting phases are shifted and favored over the p-wave solution. This fact leads to an interesting feature, namely, 
at fixed band-filling and coupling constants, decreasing the temperature from Tc defined as the critical tempera- 
ture for a paramagnet - superconductor type transition, an additional second phase transition of superconductor - 
superconductor type appears at a lower temperature < T c . This is a first-order transition between two supercon- 
ducting phases of different symmetry close to the phase boundary lines. A similar double transition between s - and 
d x 2_ y 2-wave states has been deduced by Dahm et al. in case of superconducting states emerging from spin-fluctuation 
mechanism. 

Away from half-filling (Fig. 0(d)), the increasing value of doping modifies the p-wave region of the phase diagram, 
however, the rate of paramagnetic phase emergence as the function of Tis strongly doping dependent. Furthermore, 
at a fixed value of the doping, for strong repulsive on-site interactions and not excessively low U± < values, the 
d a; 2_j / 2-symmetry phase represents not only a possible solution but it is the most stable superconducting phase among 
all the solutions of the gap equation. Qualitatively similar behaviour has been conjectured in various circumstances 
(Mierzejewski and Zielinski, 1995, Gufan, Vereshkov, Toledano, Mettout, Bouzerar and Lorman, 1995). Concerning 
the spin and charge density wave states not taken explicitly into account in Fig. [j], we expect their influence to be 
minor in the non-paramagnetic domain of the phase diagram. In this respect we would like to note that the nesting 
condition is hindered by doping in general, so its effect is reduced at low temperatures (Inaba, Matsukawa, Saitoh 
and Fukuyama, 1996) and following this line, doping usually enhances the emergence of superconductivity (Scalapino, 
1995, Iglesias, Bernhard and Gusmao, 1995). 



C. Nearest-neighbour interactions with more than nearest neighbour hopping terms 

From this paragraph on, besides NN terms, NNN contributions are also incorporated in the model Hamiltonian. 
As a first step, in order to obtain a clear image about modifications arising in the phase diagrams purely due to 
the presence of hopping terms exceeding NN distance in range, all the previously used interaction terms are kept 
unchanged but a more complex kinetic energy contribution is introduced in the present paragraph. This step also 
enables us to get an overall picture on how properties related to superconductivity are modified by DOS effects 
due to long-range hopping, without mixing these effects with the effects of NNN two-particle interaction terms. To 
accentuate the relevance to high-T c materials, we considered a dispersion relation obtained as a tight- binding fit to 
normal state ARPES data taken on Bi2212 with nonzero hopping elements up to the fifth neighbours. The numerical 
values of the used t\ terms can be found in Ref. (Fehrenbacher and Norman, 1995). Here we note, that this new, 
extended dispersion breaks the 6 — > — S symmetry of the phase diagrams and in what follows, only hole-doped cases 
are discussed. 

To show the source of modifications due to the extension in the hopping range, we compared in Fig. || the dispersion 
containing only NN, and the dispersion containing higher order contributions as well. The insets present energy 
dispersions along different directions of the first Brillouin zone. At first glance hopping terms exceeding NN distance 
in range seem to represent negligibly small corrections (for example in the present calculation |i2|/|£i| ~ 0.25 and 
|t3|/|ij.| ~0.08). This is not true, however, in the process of calculating thermodynamic averages, since their effects 
during k-integrations are extraordinarily enhanced by the supplementary extrema introduced to the dispersion. As 
it can be seen, the characteristic change occurs at small k's in the vicinity of the origin, with the flattening of the 
dispersion in this region. This leads to the fact (comparing Fig. ||(e) and (f) with Fig. ||(b) and (c)) that the region 
where JV^ 1 ~ de(k)/dk ~ holds, is much larger for the ti ^ (I < 5) case than in the ti ^ (I < 1) case. Note, 
that all the integrands within the expressions that must be evaluated are proportional to JV. As a consequence, 
main DOS effect contributions like van Hove singularities (Blumberg et al.), saddle points (Abrikosov, 1995) and 
peaks (Cappelluti and Pietronero, 1996) are strongly wiped out by taking ti = for l> 2, a choice that could lead to 
completely misleading conclusions in connection with the studied materials. 

To illustrate the effect of dispersion modifications presented in Fig. ^|. on different quantities of interest, the ratio 
R = 2Ao/ksT c , (Aq is the zero temperature gap ampitude) versus doping is plotted in Fig. |[ for a stable d x 2_ y 2- 
symmetry s uperco nducting state of the phase diagram. In this plot open triangles are calculated with the dispersion 



used in Sec. [II B| , while open circles are obtained applying the extended dispersion used in the present Section. The 



effect of hopping terms exceeding NN distance in range can easily be seen. Besides the clear absolute value differences 
we note that the curve connecting the open triangles has a minimum at about (5 = 0.1, while the other line takes its 
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least value at half filling (5 = 0.0). Here we also would like to underline that the i?~5.0 value in the optimally doped 
regime (5 « 0.12) qualitatively agrees with predictions of other authors (Abrikosov, 1995, Dagotto, Nazarenko and 
Moreo, 1995). Furthermore, Fig. |j. nicely illustrates that higher amount of doping leads to a greater value of R, an 
enhancement that is clearly increased by the ti^O (l>2) terms. 
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FIG. 2. Quaziparticle dispersion with nearest neighbour hopping terms [plot (a)] and beyond nearest neighbour hopping 
terms [plot (d)] in the first Brillouin zone. The insets show the dispersion in special directions, plots (b) and (e) are taken in 
the (-7T, 0) while plots (c) and (f) in the (7T, 7t) directions. 
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Concerning the modifications produced in the phase diagram, the obtained results are exemplified at nonzero 
temperature (ksT/W = 0.005) with Fig. [I|(a) and (b) for (5 = 0.0 and 5 = 0.2, respectively, which should be compared 
to Fig. ||(b) and (d) plotted at the same values of T and S. As it can be seen, the effects become more robust with 
increasing the doping. The extension of the Ai-symmetry domain becomes reduced and the d x 2_ y 2 phase boundary 
is withdrawn more accentuately with increasing temperature in the presence of NNN hoppings. 




0.25 



FIG. 3. Doping dependence of 2A(T = Q)/kBT c for the stable d x 2_ y 2 gap at a fixed set of coupling constants Uq = 0.2W, 
Ui = — 0.15W and XJ% = 0.0. The calculated values are denoted by triangles and circles for the cases where only nearest neighbour 
and further than nearest neighbour hopping terms were taken into consideration, respectively. 

All these theoretical results are perfectly compatible with and underline aspects observed in an increasing number 
of experimental data covering a large spectrum and suggesting the importance of NNN hopping contributions in 
building up superconducting properties: inelastic neutron scattering (Bourges, Regnault, Sidis and Vettier, 1996), 
magneto- resistance and conductivity (Mashimoto, Nakao, Kado and Koshizuka, 1996), thermopower, resistivity 
(Gasumyants, Ageev, Vladimirskaya, Smirnov, Kazanskiy and Kaydanov, 1996), and Hall effect (Hopfengartner, 
Leghissa, Kreiselmeyer, Holzapfel, Schmitt and Ischenko, 1993) measurements. Theoretical support of this line (Cap- 
pelluti et at, Normand, H. Kohno and H. Fukuyama, 1996) is also present (see also the Introduction). 



D. Influence of next-nearest neighbour interactions 



From now on, the attention is focused on the changes i n sup erconducting properties introduced by the nonzero 
value of NNN interactions. The dispersion used in Sec. |lll C| has been kept unchanged for the present Section 
and NNN coupling with U2 7^ has been considered. The obtained results regarding the superconducting phases at 
different temperature and doping values are shown in Fig. || and Fig. || for attractive and repulsive on-site inter- 
actions, respectively. Concerning the notations we would like to underline the following aspects. The presence of 
nonzero U2 term results in new singlet superconducting states characterized by k-dependent order parameters pro- 
portional to sin fcjSin k y or cosfc^cosfcy denoted by d xy or s xy , respectively, and new triplet superconducting states 
with (sin k x cos k y ± sin k y cos k x ) type of k-dependence. For the last two states the p%- and £»2-wave notations are 
used. The symmetry classification of the new states can be found in Table I. In what follows, considering the phase 
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diagrams, the notation d xy stands for an order parameter having Bi symmetry, d x ?_ y ? means an order parameter of 
Bi symmetry and the symbol A\ refers to an order parameter having s,s* and s xy components. Further notations are 
the same as in Fig. |[ In the presentation, for the sake of clarity, we concentrate on fixed values of Uq with opposite 
signs (Uq/W = —0.1 and Uq/W = 0.1 in Fig. || and Fig. ||, respectively), in order to maintain the two-dimensional 
structure of figures. 
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FIG. 4. Phase diagram at ksT — 0.005VF temperature for different dopings, using dispersion containing hopping terms 
beyond nearest neighbouring ones. Notations and curves have the same meaning as in Fig. 



The T = ground state phase diagram for different dopings (S — 0.0 and 5 = 0.1) in the presence of attractive 
on-site interaction presented in Fig. ^|(a) and (c) shows a richness of different superconducting symmetry species. At 
this point we would like to emphasize the stability of the ^i-symmetry pairing state containing a coexistence of s*, 
s xy and s-wave components. It is interesting to note that the s* and s xy components are governed by NN and NNN 
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interaction, respectively, while the conventional s-wave part is primarily determined by the on-site interaction. As a 
consequence, despite the repulsive intersite couplings, the ^-symmetry superconducting phase stabilizes due to the 
compensation of repulsions by the on-site attraction. In this compensation doping is of great importance, since a 
small increase in its value induces the emergence of Ai-symmetry pairing state in the non-superconducting region of 
the phase diagram, as shown in Fig. ||(c). 
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FIG. 5. Phase diagram with next-nearest neighbour couplings at zero (ksT — 0.0) and nonzero (ksT = 0.005W) temperatures 
for different dopings and fixed on-site attraction, [7 = — 0.1VK, for a dispersion containing also hopping terms beyond nearest 
neighbouring ones. Dotted lines mean coordinate axis and solid lines represent boundaries between different symmetry phases. 
The meaning of notations d x 2_ y 2, d xy , Ai, p-wave and paramagnetic is the same as in the text. 



The structure of the phase diagram at a higher, nonzero temperature is exemplified in Fig. ||(b) and (d). While the 
Ai phase is rapidly destroyed by increasing the temperature, the d xy phase entirely induced by NNN terms shows a 
more accentuated resistance against thermal pair-breaking effects. In this way, for large values of \U2\/W the stability 
of the day-symmetry pairing state can be coherently deduced within the presented model, a fact that underlines the 
correctness of the conjecture made by Fehrenbacher and Norman (Fehrenbacher et al. related to the stability of the 

dxy st&tc 

Concerning triplet pairing we note that the states with p-, p2- and ^-symmetry never have lower free-energy in the 
coupling constants' space under study than the triplet pairing state of standard p-wave symmetry. 

If now the Uq > 0.0 case is considered, the phase diagram significantly changes, as it is shown in Fig. ^. The most 
important feature of this figure is, that the Ai-symmetry solution never reaches the level of a stable superconducting 
state within the analyzed U2 region for |J7i|<Cl.O for slight dopings. This result is in accordance with the observation 
of Fehrenbacher et ai, where it is claimed that in the presence of on-site repulsion with small NNN attractions 
(\U2\/W <g. 1.0), the stable superconducting state of Ai-symmetry can be realized in the limit of |i7i| — > 0.0 only at 
relatively high values of the doping. 

Figure ||(a) and (c) show the T = phase diagram of the model with Uq = 0.1W at two doping values ((5 = 0.0 and 
5 = 0.1) in the (Ui/W, U2/W) plane. The phase boundary between the d x 2_ y 2 and d xy phases moves slowly into the 
domain of d x 2_ y 2 phase with increasing values of the doping, that is, higher doping favors states of d^y-symmetry over 
states of c? a .2_ y 2-symmetry. Furthermore, as it can be observed in Fig. ^(a) and (c), the d x 2_ y 2-d xy phase boundary 
is linear and strongly affected by doping. For small values of 5 it coincides with the Ui =aU2 line with a ~ 1.0, while 
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for larger doping values it is shifted away and a becomes smaller than 1.0. This leads to the following interesting 
feature: changing the doping continuously the symmetry of the stable superconducting state changes and a transition 
between the two different d-like pairing states (from £?2-symmetry to Bi-symmetry) takes place for small, fixed values 
of the coupling constants. This transition line within the T = phase diagram, suggested by detailed investigations, 
is of first order. 
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FIG. 6. Phase diagram with next-nearest neighbour couplings at zero (fcgT = 0.0) and nonzero (ksT = 0.005W) temperatures 
for different dopings and fixed on-site repulsion, Uo = 0.1W , for a dispersion containing also hopping terms beyond nearest 
neighbouring ones. Notations and curves have the same meaning as in Fig. n. 



Examining phase boundaries between the p and d xy or d x 2_ y 2 superconducting states, one can see them moving 
upward with increasing doping into the d a ,2_ y 2-symmetry phase. However, large NNN interaction (|£/2|/W^ 1) and 
small AW attraction (|fi|/W^C 1) favors superconducting pairing state of d x i_ y i or d xy symmetry, depending on the 
sign of U2, an observation that holds even for heavily doped systems. 

Turning now to the T ^ case depicted in Fig. ^(b) and (d) , it can be observed that the triple-point moves downward 
in the U\, U2 < 0.0 region with increasing the temperature. The phase boundary of the d x 2_ y i phase modifies more 
rapidly with increasing T, suggesting that in the large | tT" 2 1 region (with U2/W ^ < 0.0 and \U\\/W 1.0) the highest 
critical temperature can be reached with a d xy phase. The stability of the d xy phase against thermal fluctuations, in 
comparison with the same stability of d x 2_ y 2 phase, is again striking. 

In the light of Figs. ^| and ^, one can conclude that small |f7i|/W and large attractive NNN coupling U2 favors d xy 
pairing, while small IC/2I/VF and strong attractive AW coupling U\ gives rise to stable d x 2_ y 2 pairing. These results 
emphasize the correctness of the conjecture made by Wenger and Ostlund (Wenger et al.) related to the stability 
of different d symmetry species, who also claimed that third-neighbour attractive interaction could lead to a stable 
k-dependent s-wave within the phase diagram, especially in the presence of repulsive closer-neighbour interactions. 

The electron doping versus hole doping behaviour is examplified with Fig. [?], where doping dependence of gap 
amplitudes and free-energy are shown for a certain set of coupling constants. Having a glance at Fig. 0(a) one can see 
that besides the gap amplitudes, phase also changes (i.e. sign change of Aj in the presented figure) emerge for different 
solutions at a critical value of doping. More interestingly, Fig. 0(b) predicts the stable superconducting state to have 
Ai-symmetry in the electron doped regime even for repulsive on-site couplings. This suggests that the symmetry of 
the pairing state in electron doped materials (for example in NCCO) might be explained without assuming attracion 
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in the on-site interaction channel. Systematic study of the electron doped case is under progress. 
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FIG. 7. Doping dependence of gap amplitudes [plot (a)] and of corresponding free-energies [plot (b)] at a fixed set of coupling 
constants Uo = 0.1W , Ui = — 0.1W, XJi =— 0.2W^ for the dispersion containing hopping terms beyond nearest neighbouring ones. 



In Fig. U we focus our attention on the effects accentuately caused by NNN interaction terms in T c , in zero 
temperature gap-amplitude and in the ratio of R = 2Ao/ksT c . For simplicity, U\ is set to be zero in this study. 
This particular choice results in an order parameter of the form A(fc) = A s -(-A Sx cos fc^cos k y , (i.e. only s and s xy 
components). The values of Uq and U2 are tuned so that the emerging Ai-symmetry pairing state becomes the stable 
superconducting state for the chosen set of interaction parameters. The effect of U2 on the mentioned quantities is 
presented in the figure. We note, that for Ui ^ 0.0 values a similar analysis can be done leading essentially to the 
same features. 

First of all, Fig. ||(a) shows the t/ 2 dependence of the gap amplitudes at zero temperatures. One can immediately 
realize the exponential type behaviour of A Sxb versus — U2 for small, IL^I/PF^ 1.0 values. At the same time, A s 
has a linear dependence on — £/ 2 . The U2 dependence of the critical temperature is exemplified in Fig. ||(b) for the 
Ai-symmetry solution. A good fit to the numerically calculated points can be obtained by the T c ~ exp[— l/(Ky/U2~)] 
relation, where K = 0.5226. Such type of T c behaviour related to non-on-site interactions is not unusual. For example, 
a similar relation was suggested by Dahm et al. (Dahm et al.) for the nearest neighbour coupling dependence of T c , 
instead of an exp[— \/{Kgi)\ functional form. Concerning the ratio of the zero temperature gap amplitude to T c , in 
Fig. ||(c) the R — 2A /kBT c values are plotted versus — U2/W for the two components of the Ai solution. In the 1^ = 
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case (where the order parameter has the form of A(fc) = A s ) this ratio is about 3.52 and hence corresponds to the 
BCS predicted value. However, if NNN interaction is turned on an s xy component of the order parameter develops 
rapidly resulting in an increase (decrease) of the ratio R for the s xy (s) channel. 
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Based on these results it can be seen that nonzero next-nearest neighbour coupling, besides the emergence of 
different possible superconducting species, also influences the main characteristic properties of the superconducting 
phase. 

IV. SUMMARY AND CONCLUSIONS 

The superconducting phase diagram of the two-dimensional extended Hubbard model containing hopping and 
interaction terms with spatial dependences exceeding nearest neighbour distance in range has been systematically 
analyzed in mean-field approximation. The possible superconducting phases of different symmetry were taken into 
account unrestrictedly. For building up the phase diagrams, the emerging phase was choosen on the basis of a free- 
energy (or at T = ground state energy) analysis of all the possible solutions of the gap equations for every fixed 
set of coupling constants, temperature and doping values. Main characteristics of different superconducting phases 
including critical temperatures T c , zero temperature gap amplitudes Ao, Aq/T c ratios, temperature and doping 
dependences were also studied in detail for every domain of the phase diagram. The obtained results clearly underline 
the importance of next nearest neighbour terms (both hopping and interaction contributions) in determining the 
prominent superconducting properties of the system. 

The present paper explored such a parameter domain of the phase diagram which has not been investigated so far 
systematically in literature starting with Fig. where in the presence of only nearest neighbour terms, non half-filled 
band case was also investigated. Starting from this level different next-nearest neighbour contributions were taken 
gradually into consideration within the extended Hubbard Hamiltonian, their effects being systematically presented 
in Figs. |^-^ and analyzed in detail with physical implications being emphasized, especially those connected to the 
presence of next-nearest neighbour terms. 

On this line the following interesting features deserve attention: double phase transitions with decreasing tem- 
perature (following the emergence of a superconducting state at T c the symmetry of the superconducting order 

parameter changes at rj 2 ^ < T c in a first-order phase transition); quantum phase transitions at zero temperature 
between pairing states of different symmetry driven by doping; enhancement in the resistance of the stable supercon- 
ducting phase against thermal pair-breaking effects in the presence of interaction terms exceeding nearest neighbour 
distance in range; increase in the T = gap-amplitude over T c ratio for some symmetry species; the main supercon- 
ducting properties become more sensible with respect to doping; strong asymmetry between hole and electron doped 
cases; rich spectrum of stable superconducting states of different symmetry (in the presence of repulsive on-site and 
in some cases even for repulsive nearest neighbour interactions, the A\ symmetry is favored for electron, and a d-wave 
type order parameter for hole doping values); elimination of k-independent order parameters even at half- filling, and 
so on. 

We strongly hope, that the presented results will constitute a valuable starting point for further theoretical investi- 
gations related to the effect of long-range contributions in building up superconducting properties of the system under 
study. 
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